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Abstract 

Following an earlier derivation by Catani-de Florian-Grazzini (2000) on the scheme dependence in the Collins-Soper- 
Sterman (CSS) resummation formalism in hard scattering processes, we investigate the scheme dependence of the 
Transverse Momentum Distributions (TMDs) and their applications. By adopting a universal C-coefficient function 
associated with the integrated parton distributions, the difference between various TMD schemes can be attributed to 
a perturbative calculable function depending on the hard momentum scale. We further apply several TMD schemes to 
the Drell-Yan process of lepton pair production in hadronic collisions, and find that the constrained non-perturbative 
form factors in different schemes are remarkably consistent with each other and with that of the standard CSS formal¬ 
ism. 


1. Introduction 

The Transverse Momentum Distributions (TMDs) and the nucleon tomography in momentum space have attracted 
strong interest in recent years iiiiia. TMDs provide a unique opportunity to investigate the novel correlations between 
the parton momentum and the nucleon spin. They unveil the strong interaction QCD dynamics in a manifest way, such 
as the gauge invariance leading to the sign change Ejia of certain TMDs in different hard scattering processes, and 
the QCD factorization and evolution which are crucial for predicting the scale dependence of the spin asymmetries. 
On the theory side, the TMDs are not straightforward extensions la of the conventional collinear parton distribution 
functions (PDFs). They hold special properties that differ them from collinear PDFs and play important roles in high 
energy scattering. The associated phenomena are direct consequences of perturbation gauge theory computation of 
the famous Sudakov form factors lb) back in 1950s. 

When one applies the TMD factorization to physical processes, one has to consider the associated QCD dynamics 
in the definition of TMDs and in the phenomenological studies. Especially, large logarithmic corrections from high 
order perturbative calculations have to be taken into account and resummed Q to all orders. In addition, the naive 
gauge invariant TMD definition contains the so-called light-cone singularities at higher orders and needs to be regu¬ 
lated ||8l. Several ways to implement such a regularization have been proposed in the literature and they introduce the 
scheme dependence in TMDs and their applications ||8]|9][T0][TTl[T2l|T3l. The goal of this paper is to investigate such 
scheme dependence, which is of crucial importance for applying the TMDs in hard scattering processes and extracting 
the associated nucleon structure from experiments. 

In the context of the standard Collins-Soper-Sterman (CSS) resummation formalism 13 , the TMDs are expressed 
in terms of the collinear parton distributions via an additional factorization at small b <Si 1/Aqcd, where b represents 
the Fourier conjugate variable associated with the transverse momentum kj_. The final expressions for the measured 
cross-sections differential in transverse momentum of the observed particles do not depend how we define the TMDs 
at the first place when such relations to collinear PDFs are used. In other words, in all TMD formalisms of Refs. El 
ismoiiiiiiiiiini, one will obtain the same results as that of the standard CSS resummation. However, as discussed 
in an early paper by Catani, de Florian, and Grazzini lfT4l . even in conventional CSS formalism there is freedom to 
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separate the so-called hard factor, H, which depends on the running coupling at the hard momentum scale Q from the 
C-coefhcient functions associated with the integrated parton distribution functions where running coupling depends 
on Hb = co/b with co = It was referred to in Ref. lfT4ll as the scheme dependence of CSS resummation. The 

relation between different schemes was further demonstrated by an order by order proof lfT4l . The relevant derivations 
with explicit results up to next-to-next-leading order for Drell-Yan, Higgs boson, di-photon production processes have 
been extensively discussed in Ref. ifTSlI . The same argument applies to the scheme dependence in the TMD formalism 
as well E1121 mill] [El [El. By adopting a universal C-coefhcient function associated with the collinear parton 
distributions na, the connections between diherent schemes can be attributed to the hard coefficients and can be 
established order by order in perturbation theory. As a result, all the TMD scheme dependence can be accounted for 
and the schemes can be unihed and compared to the standard CSS resummation in description of the experimental 
data in phenomenological studies. 

Furthermore, this unihcation provides an attractive interpretation for the CSS resummation, from which we have a 
clear TMD interpretation of the hard scattering processes. To establish this, we apply this scheme in the global analysis 
of the Drell-Yan process of lepton pair production inpp collisions, and fit the associated non-perturbative form factors. 
In the calculations, we adopt the so-called h*-prescription and derive the relevant perturbative coefhcients following 
the procedure of Ref. IT4l . Our results show that the non-perturbative form factors are remarkably consistent with 
that in the standard CSS scheme. This is a very important result that will clarify the confusing issues of the TMD 
evolution which has been intensively explored in recent years. 

The rest of the paper is organized as follows. In Sec.we briefly introduce the TMD schemes in hard scattering 
processes, and derive the relevant coefficients. In Sec. 1^ we ht the experimental data of Drell-Yan type of hard 
processes in hadronic collisions and constrain the associated non-perturbative form factors. And finally, we conclude 
our paper in Sec. 


2. TMD Schemes 


Let us start with the standard CSS resummation formalism for Drell-Yan lepton pair production processes at low 
transverse momentum: A{Pa) + B{Pb) —> J*iq) + X ^ {'^ + (~ + X, where Pa and Pg represent the momenta of the 
incoming hadrons A and B, respectively. The differential cross section can be written as 171 . 
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where q^^ and y are transverse momentum and rapidity of the lepton pair, respectively, = Analj^jiNcsQ^ with 

s - {Pa + Pb)^- In the above equation, the first term is dominant in the q_i_ Q region and Wuu denotes the all-order 
resummation result which has the following form QHl: 
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where pb - co/h* with co = 2e“^^ and the Euler constant, xi _2 = I represent the momentum fractions 
carried by the incoming quark and antiquark in the Drell-Yan processes, the symbol ® for convolution in xi(x 2 ) and 
flA(x,pb) and fj/Bix,pb) stand for the collinear integrated parton distribution functions at the scale pb- In Eq- 0, 
-prescription, b ^ b, - bj is introduced m. The form factor S {Q, b) contains perturbative and 

nonperturbative parts, such that the total form factor for quarks can be written as S (Q, b) - SperiiQ, bP) + Snp{Q, b). 
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where A, B and C coefficients calculable order by order in perturbation theory perturbative series A - (aslTT)", 

B = (asln)”, C - (asln)”. The A, B, C coefficients can be derived llT4l . 
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in the standard CSS scheme. In the standard CSS formalism, the hard coefficient HcssiosiQ)) = 1 for all orders. 

We would like to emphasize that the resummation formula and the associated coefficients are uniquely determined, 
once the scheme is fixed. The reason is simple. In the perturbative calculations of hard processes at low transverse 
momentum, the large logarithms depend on two separate scales: Q and l/b, the hard momentum and the Fourier 
conjugate of the traverse momentum q^, respectively. The resummation of these large logarithms has to take the 
form as in Eq. 0, as a consequence of perturbation gauge theory computation of Sudakov form factors ElIllEl. 
Additional factors in the CSS resummation come from the fact that the collinear gluon splitting is proportional to 
l/q]_ (again a result of a gauge theory computation), for which the Fourier transformation leads to ln(//*/yu) where /i 
represents the PDF scale. Therefore, the integrated parton distribution is calculated at for canonical choice of the 
resummation. By doing so, we also resum the logarithms associated with collinear gluon radiation. The coefficients A, 
B, and C can be obtained from the factorization derivation, or by comparing to the fixed order perturbative calculations. 
For phenomenological applications, the CSS formalism has been very successful in Drell-Yan lepton pair production, 
W^/Z boson production in hadron collisions lfni[T8lfT9ll20ll2T1l22lE3ll24ll^l^l27ll28ll29ll. 

As discussed in Ref. 03, there is a freedom to absorb as(Q) corrections from higher orders in the definition of 
hard coefficient H^^^\asiQ)) of Eq- §. Then, the associated B and C coefficients will be modified according to the 
renormalization group equations. This was referred to as the scheme dependence in the CSS resummation in Ref. 03- 
In the following, we will apply this idea to discuss the TMD interpretation of the CSS resummation formalism, where 
the scheme dependence is essential in the TMD definition and factorization. 

TMD factorization ||9] aims at separating well defined TMD distributions in Eq. Q, such that the TMD distribu¬ 
tions can be used in different processes in a universal manner. The fe-space expression of Wuu in Eq. 0 thus can be 
rewritten in the TMD factorization in terms of a product of process independent TMDs and a process dependent hard 
factor: 

where both the subtracted quark distribution and hard factor depend on the scheme we choose to regulate 

the light-cone singularity in the TMD definition. In tfiis paper we consider three TMD schemes: (1) Ji-Ma-Yuan 
2004 mnillll; (2) Collins 2011 0; (3) Lattice IHl or Collins-Soper 1981 0. The so-called FIS scheme was 
shown to be equivalent to Collins 2011 scheme ll30l . Moreover, because of usage of space-like gauge fink in the 
lattice scheme, the results in this scheme coincide with the the original Collins-Soper 81 scheme. Extensions to other 
formalisms can follow accordingly. 

We take the example of Ji-Ma-Yuan 2004 scheme (JMY) IfTOl [TTII . where the unpolarized quark distribution is 
defined as 

with the gauge fink Xi.(-°°;f) = exp(-/g^ dAv ■ A{Av + The above definition contains the light-cone sin¬ 
gularity if we take the gauge fink along the light-front direction, = 0. The way to regulate this singularity and 
subtract soft gluon contribution defines the scheme for TMDs. In the JMY scheme, the gauge fink is chosen to be 
slightly off-light-cone, such that n = (l“,0^,0x) v = (v“, v^,0x) with v~ » Similarly, for the TMD antiquark 
distribution, v was introduced, v = (v', v^,0j_) with v* » v~. Because of the additional v and v, there are additional 
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invariants: = (2v ■ Pa?'I v^, - (2v ■ PeYI v^, and - (2v ■ v)^/v^v^. Accordingly, the soft factor is defined as, 

5 "’Hb) = <oi 4(Z7 jxixfe ji:„(0)x,(0)io). (?) 

Following the subtraction procedure of Ref. a, we can define the subtracted TMDs in Z7-space in the JMY scheme 
as. 


Asub) 
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where fq{x,b',(,pF,p) is the /j-space expression for the un-subtracted TMD of Eq. The evolution equations are 
derived for the TMDs: one is the energy evolution equation respect to the so-called Collins-Soper evolution equa¬ 
tion 0 and the renormalization group equation associated with the factorization scale pf related to anomalous 
dimensions of the distribution /. After solving the evolution equations and expressing the TMDs in terms of the 
integrated parton distributions to have a complete resummation results, we can write. 


= Q,P) 


®fi(x,ph), 
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where we have chosen the energy parameter = pQ^ and the factorization scale pF - G to resum large loga¬ 
rithms nniin]. The perturbative form factor 5^^^, contains contributions from the Collins-Soper evolution kernel 
and the renormalization equation respect to the factorization scale pF and ph- Similar to the CSS resummation, /?»- 
prescription was applied. In the above equation, we have also followed the derivations of Ref. flTl to include the 
p-dependence in the hard factor by applying the renormalization group equation of running coupling ccj. By doing 
that, the C-coefficients are much simplified and have the following universal TMD form[^ 
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for the quark-quark and quark-gluon splitting case. A universal C-function in the CSS resummation formalism has 
also been emphasized in Ref. ca. From the results in Ref. IfTOl . see, for example, Eq. (36) of llTOl . we obtain 


(a,(G);p) = l + ^Cf(lnp 


In^p 



( 12 ) 


The above equations are derived based on the perturbative calculation and the associated QCD factorization for the 
TMDs. They apply to all TMD schemes ISJ 0 [TOl [TT] [121 El mentioned above. The collinear divergence in the 
TMDs can be factorized into the integrated parton distributions as shown in Eq. 0 at small b «c l/Agco- For large 
b, a non-perturbative function has to be included. The universal C-coefficient function is adopted to simplify the final 
expression for the TMDs and minimize the higher order cotTections associated with the integrated parton distributions. 

Similarly, for the Collins 2011 (ICC) scheme, we have ll9l f3T1l . 






i 

(13) 


= 1+0(4). 
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*In principle, we can also choose Cess of Eq. for the C-coefficients, which will go hack to the standard CSS resummation for phenomeno¬ 
logical applications. We chose these coefficients for simplicity. 
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where C,c is the regulation parameter in JCC scheme and the as correction in vanishes ^ Again, we emphasize 
that C-coefficient takes the same form as that in Eq. ( [T0] i. Therefore, the scheme dependence in the TMDs only comes 
from the hard function Tq as we have shown in the above equation. 

Recently, there has been a motivated study to formulate the TMDs on lattice, where a different subtraction scheme 
was adopted, for which we have ini 

4(La!.M’ = e", /UF = 0) = ^Lat. 

(15) 

T^^'iasiQ)) = l + gCf(-2), (16) 

where the regulator f is defined as = (2«j • P)^l(-n^) with space-like n^: - -1, n^ - P - -P^. As we mentioned 

above, lattice scheme uses the same space-like gauge link as the original Collins-Soper 1981 scheme, that is why the 
final expression for TMD are the same in both schemes. 

Applying the above TMDs into the factorization formula of Eq. 0, and comparing to that in Eq. we will 
find that the TMDs actually provide a special scheme for the CSS resummation in the context of Ref jUj. We can 
immediately derive the relevant coefficients, 

4md («^(0)) = («^(0)) X («.^(0)) X 'Hjm) (0; 0) - (17) 

which will enter into Eq. (El) for phenomenological applications. Because the C-coefficients are universal among 

•——I (DY^ 

different TMD schemes, we conclude that will be the same in all of the three schemes discussed above. In 

particular, in the JMY scheme, all three factors in Eq. •EH depend on p, however the final result for does not 

depend on p. This demonstrates that all the TMD factorization schemes are equivalent in the context of the CSS 
resummation formalism, which will be used in the phenomenological applications. This can be verified from the 
above explicit results and from the associated hard factors calculated for different schemes at the one-loop order, and 
order by order proof can be done accordingly. 

Eurther comparison also indicates that the perturbative and non-perturbative form factors for the quark and anti¬ 
quark can be related to that in the CSS formalism Eq. 

S%{Q, K) - K) = SperpiQ, K)I2 , (18) 

SUQ^b) + Sl^(Q,b)^S^p(Q,b), (19) 

where the perturbative form factor S pert(Q,b,) takes the form of Eq. (j^ with A and B coefficients for a particular 
TMD scheme. The above equation for the perturbative form factors can be verified explicitly from one-loop results in 
the TMD factorization of Refs. Il9l ITOl [TTl fT3l . Higher orders can be calculated in perturbative expansion. 

Erom the above one-loop results for 'Fq,^ and the relevant hard factors in the TMD factorization calculated in 

Refs, nniiniisiiii, 

^^^(2;p) = 1 + ^Cf (3 In ^ - In" ^ + tt" - 8j , (20) 

,^™djmY( 0;^) ^ J +Inp2)in^ -ln2pH-ln"pH-2;r" -4j , (21) 

■K™DLat.(g;^) ^ j ^ ^Cf (in ^ + tt" - 4j , (22) 

we obtain the one-loop expression for //tmd as, 

77™^ = \cF{n^-^), (23) 


^There is an ambiguity for the ultra-violet (UV) subtraction: an additional term of n-fll should be added in as correction if we follow the 
standai'd MS subtraction used in the standard CSS. Here we adopt Collins-11 prescription for the UV subtraction. 
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so that 


4MD(e) = i + 

For B and C coefficients, following the derivation of Ref. m, we will obtain 
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whereySo = A*'^) and B*'* remain the same as the standard CSS scheme. We will apply these coefficients 

in the next section to analyze the Drell-Yan type of lepton pair production in hadronic processes to constrain the 
associated non-perturbative form factors. 

Following the arguments of Ref. lfT4l . the process-dependence is included in H in Eq. (|^, so that C-coefficients 
will be universal. We can apply the same C-functions to the quark distributions in other processes, such as the Semi- 
Inclusive Deep Inelastic Scattering (SIDIS), for which we have the standard CSS resummation coefficient 

5 - /3o Cf 7? 12 , (26) 

with front Eq. fl and all other coefficients that have been listed in Refs, ll^ |28]| . The hard function for 
SIDIS is 


H. 


{SIDIS) 

TMD 


asiO) 

(0 = l + ^pCf(-8) 

2n 
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if we choose the TMD scheme for this process. 


3. Non-perturbative Form Factors and TMD Interpretation 

As we mentioned above, we will apply the f>t-prescription for the non-perturbative form factors. We will follow 
the SIYY parameterization ||29]| . This parameterization is motivated by a phenomenological study Il28l and is inspired 
by matching to perturbative calculations of the Sudakov form factors 01 Ea. It has the following form, 

5np(0^) = g2\HblK)\n{QIQo) + gib\ (28) 

with the initial scale - 2.4 GeV^ and cut-off parameter b„jax =1.5 GeV”'. The parameters gi 2 have been fitted 
to the experimental data of Drell-Yan type processes in Ref. Il29l using the standard CSS formalism (H = 1). In the 
study of Ref. Il29l it was found that the experimental data are consistent with jr-independent non perturbative factors. 
In the following studies, we will take the above simple form of Eq. ( |28] l. Since we will compare the TMD schemes to 
the standard CSS scheme, we will keep all relevant parameters in non-perturbative factors fixed except for the changes 
in the coefficients H‘'^\ and 

We compare our results to the same set of the experimental data sets as those used in Ref. Il29l . The data sets 
include the Drell-Yan lepton pair production from fixed target hadronic collisions from R209, E288 and E605 ll3^ 
|37l|38l, and Z boson production in hadronic collisions from Tevatron Run I and Run II ll3^ l40l HTI |4^ . We proceed 
with the fit of the experimental data using standard CSS, done in Ref. Il29ll and the TMD-schemes described in this 
paper. Notice that all TMD-schemes have exactly the same hard factor 7/^^, so by doing a single fit we effectively 
obtain underlying TMDs in either Collins 2011 ID, Lattice lfT3l . or Ji-Ma-Yuan 2004 IfTOl [TTI schemes. The fitted 
parameters are found to be. 


SIYY®: 

SIYYtmd: 


0.212, g 2 = 0.84 , total ^ 

0.212, g 2 = 0.84 , total = 168 , 


(29) 

(30) 


where the first line is for the standard CSS scheme fit 11291 . the second for the TMD-scheme JMY and ICC with 
coefficients in Eqs. ( 10|11|23|25 1. There is hardly any difference between the two fits. There is no difference in the 
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Figure 1: TMD up-quark distributions = 0.1,/?) as functions of b at different scale = 2.4, 10, 90 (GeV^) for three different schemes, 

from the top to the bottom: JCC (§), Lattice (H), JMY Gnillll (Inp = 1). 
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Figure 2: TMD up-quark distributions = 0.l,k±) as functions of the transverse momentum kj_ (GeV) at three different scales = 2.4, 

10, 90 (GeV^) for three different schemes, from the top to the bottom: JCC (9l, Lattice □3, JMY (T0][TT1 (lnp= 1). 
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comparisons to the experimental data either. It demonstrates the effective equivalence between all schemes in the 
phenomenological studies. Theoretically, the difference could come from higher orders, such as a] in //tmd and the 


coefficients at N^LL for the resummation which are beyond what we have considered in this paper and Ref. 

With the non-perturbative form factors determined from the experimental data, we can compare the TMD quark 
distributions in different schemes by evaluating them using Eqs. 1^ Ji-Ma-Yuan 2004 scheme has a residual 
dependence on the value of p, we fix it by choosing Inp = 1 .The transverse momentum dependence in three schemes 
is calculated by Fourier transformation respect to b using Eqs. ( 9|13|15| l. In Fig.[T] as an example, we plot the up-quark 
distributions = 0.1,/?) at .r = 0.1 for different schemes at different scale - 2.4, 10, 90 (GeV^) as functions 

of b (GeV ') and in Fig.|^we plot - 0.1, as function of the transverse momentum One can see from 

Fig .[^that at low values of the non-perturbative part of the distribution becomes very important and the values of 
b > bmax dominate the result in space. At higher values of the large b tail of the distribution is suppressed and 
the whole distribution can be computed using mainly perturbative regime b < b^ax- In this regime the results will have 
a relatively low sensitivity to the non-pertu rbative inp ut of TMD evolution. Again, the difference between different 


schemes is due to the coefficient 'Fq in Eqs. (12 


14 


16 1 . Because the difference is proportional to as(Q), it will become 


smaller at higher scale Q as shown in Figs.[]J|2 Similar plots have been shown in Ref. OTI for the quark distributions 
in the JCC scheme, however, using the previous BLNY parameterization uni for the non-perturbative form factors 
obtained in CSS resummation. In our calculation we have consistently used relation between different schemes and 
the fitted non-perturbative form factors. 


4. Conclusion 

In this paper, we have investigated the scheme dependence in the TMD paiton distributions and factorizations to 
describe the experimental data of hard scattering processes in hadron collisions. The equivalence between different 
schemes can be proven in perturbation theory order by order following the procedure of a similar study of Catani- 
de Florian-Grazzini 2000 03. We have studied three such schemes, Collins 2011 a, Lattice ISl, or Ji-Ma-Yuan 
2004 EOlini, and have demonstrated the equivalence between them and equivalence to the standard CSS method. 
The associated coefficients are illustrated at one-loop order. 

With TMD scheme dependence embedded in the coefficients, F and the resummation formulas have been 

applied to the Drell-Yan type of lepton pair production in pp collisions, and the associated non-perturbative form 
factors are determined from the global fit. We found that the TMD-schemes produce the same phenomenological 
results as compared to the standard CSS scheme for the resummation. More importantly, the parameters of the 
associated non-perturbative form factors are also found to be the same in all schemes. 

Using the fitted parameters, we can calculate the TMD quark distributions as functions of the transverse momen¬ 
tum. We have compared the results from three different schemes. This comparison becomes useful in the TMD 
interpretation of the experimental results. 

In this paper we explored the spin-average quark distributions. Similar studies can be carried out for all other 
TMDs. In particular, the quark Sivers function, which describes the correlation of the transverse momentum of the 
quark and the nucleon spin, can be formulated in the CSS resummation formalism. The non-perturbative form factors 
for Sivers function, however, will be different from the unpolarized quark distributions discussed in this paper. In 
order to determine these non-perturbative factors one needs to perform a global fit to the existing SIDIS data. We 
leave that for a future publication. 
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